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Abstract 

Based on the results of accurate variational calculations we demonstrate stability of the five-body 
negatively charged ions a^6"*'/x~e~e~. Each of these five-body ions contains two electrons e~ , one 
negatively charged muon /i~ and two nuclei of the hydrogen isotopes a,b = (p,d,t). The bound 
state properties of these five-body ions, including their hyperfine structure, are briefly discussed. 
We also investigate the hyperfine structure of the ground states of the four-body muonic quasi- 
atoms a~^b^ fi e . In particular, we determine the hyperfine structure splittings for the ground 
Q-f state of the four-body muonic quasi-atoms: p^d^fi^e^ and p^t^fj,~e~ 
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I. INTRODUCTION 

In our earlier work [l| we have discussed the four-body muonic quasi-atoms a^b^ij,~e~ 
and some of their properties. Here and everywhere below, the notations a and b stand for 
the nuclei of the hydrogen isoptopes (p designates protium, d is deuterium and t is tritium), 
fi^ (or /i) and e~ (or e) designate the negatively charged muon and electron, respectively. 
These muonic quasi-atoms are formed in liquid hydrogen/deuterium during muon- catalyzed 
fusion. It appears that the four-body muonic quasi-atoms, e.g, p'^d'^ fi~e~ , d'^d'^ fi~e~ , etc, 
have many extraordinary properties which follow from their two-shell cluster structure. In 
particular, each of these ions has a heavy central nucleus abfi which is positively charged. 
Spatial radius of such a three-body nucleus is ~ One electron moves around this 
heavy nucleus at disctances i?e ~ cto, where gq is the Bohr radius. It is clear that R^, 
since ^ ^ 206.768. 

me 

The two-shell cluster structure of these quasi-atoms allows one to predict a large number 
of their bound state properties without extensive numerical calculations. To check the overall 
accuracy of such predictions in |l| we have performed variational computations of the ground 
states in the six muonic atoms: pp/ie, pci/ie, pt/ie, dd^e, dtfie and ttfie. In such calculations 
we used the following particle masses nip = 1836.152701 me, = 3670.483014 me, m^ = 
5496.92158 me and = 206.768262 me. These values of the particle masses are often used 
in recent computations of muonic molecular ions abfi and other related systems. The same 
masses will be used in this study. 

In this study we continue our analysis of the four-body quasi-atoms a'^b'^fi~e~ and closely 
related five-body negatively charged ions a~^b^ fi^e^e^ (or a^b^ fi^ e^)- As follows from 
our preliminary evaluations (see below) the correspoding five-body negatively charged ions 
a'^b^ fi^e~e~ are stable and can be observed in actual experiments. These ions are of great 
interest for general theory of bound states in few-body systems with unit charges. In this 
study we want to investigate the stability of these five-body ions formed by the particles 
with unit charges. Another aim is to determine the hyperfine structure splitting of the 
ground states in the five-body abfie2 ions and four-body abfie hydrogen-like quasi-atoms. 
In particular, we discuss the hyperfine structure and determine the hyperfine structure 
splitting for the ground state of the pdfie quasi-atom and predict the results of possible 
future experiements. 
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II. BOUND STATES IN THE FOUR-BODY MUONIC QUASI- ATOMS 

In our variational calculations of the four-body quasi-atoms a'^b'^j2~e~ in [l| we have 
applied the variational expansion written in the basis of the six- dimensional (or four-body) 
gaussoids of relative (or inter-particle) scalar coordinates Vij. This variational expansion 
was originally proposed thirty years ago in [s] for accurate variational calculations of few- 
nucleon nuclei. For the bound S{L = 0)— states the variational anzatz of six-dimensional 
(or four-body) gaussoids in relative coordinates is written in the form 

1 ^ 

^L=o = ^(1 + 5abtabVab) '^k " exp{-aSrl^ - aigV^g - afjrl^ - a~Srl^ - aSrl^ - af^rl^ll) 

where Ck are the linear coefficients (or linear variational parameters), while af^'' are the 
optimized non-linear parameters. The notation eabVab means the appropriate symmetrizer 
(or antisymmetrizer) , i.e. a projection operator which produces the wave function with the 
correct permutation symmetry in those cases when a = b. This case is designated in Eq.(|3]) 
with the use of the delta-function. The operator Vab is the pair-permutation operator for 
all coordinates, i.e. for the spatial, spin, iso-spin, etc, coordinates. In these calculations 
we have used two trial wave functions with = 400 and = 600 terms, respectively. All 
non-linear parameters a^j^ in Eq.(|3]) have been varied carefully. 

The results of new variational computations of the ground states in the six muonic quasi- 
atoms ppfie, pdfie, ptfie, ddfie, dtfie and ttfie can be found in Table I, where the total energies 
of these quasi-atoms are given in atomic units h = l,me = l,e = 1. These energies are 
slightly better than the analogous total energies obtained in In this study we need the 
improved total energies of the four-body muonic quasi-atoms a"'"6+/i~e~ to show the stability 
of the corresponding five-body ions a~^b'^fi~e~e~ which contain two electrons (see below). 
Note also that there are twenty two bound states in the six three-body muonic molecular 
ions pp^,pdfi,ptfi,ddfi,dtfi and ttfi. Formally, each of these ions (in their ground and/or 
excited states) can bind an electron and form the four-body quasi-atom a6/xe. In this study 
we restrict ourselves to the analysis of the six ground states only. 



3 



III. STABILITY OF THE FIVE-BODY MUONIC IONS 



As follows from the general theory, the stabihty of the four-body quasi-atoms a'^b'^fj,~e~ 
means that the corresponding two-electron ions a'^b'^ /i^e'e" (or a^b'^ ^'e^) are also stable. 
Each of these ions must have only one bound state which is, in fact, a weakly bound state, i.e. 
its binding energy is < 1 % of its total energy. Here, the stability means the energy stability 
against the following two-particle dissociation into quasi-atom abfie and free electron 

a^b'^ li~e~e~ — a^b'^ ii~e~ + e~ (2) 

Briefly, the word 'stabihty' means that the total energy of the five-body muonic ion 
a^b'^ li~ e~ e~ is lower than the total energy of the four-body a^b^fi^e^ quasi-atom in its 
ground state. Furthermore, the general theory of bound states in few-body systems predicts 
that the five-body a^b'^ ^~ ion has the only one bound S{L = 0)— state and this state is 
a singlet electron state. Formally, such a state can be designated as the state, but in 

some cases this notation can lead to a confusion. In this Section we want to prove that the 
five-body negatively charged muonic ions a'^b'^ ii''e~e~ are stable by performing accurate 
variational computations. 

First, note that the bound (ground) state of the a^b^ fi~e~e~ ion has the two-shell cluster 
structure. The internal three-body system, i.e. nucleus with electric charge +1, is formed by 
the three heavy particles a, b and /i. The spatial radius of such a 'quasi-nucleus' is i?^ ~ 
where is the muon mass and oq is the Bohr radius. Everywhere below in this study 
we shall use only atomic units. In these units ^ = l,e = 1 and mg = 1 and one finds 
= 206.768262 and ao = 1. This means that i?^ fa 206 768262 ^ '^o- electrons are 

moving around the central, positively charged 'quasi-nucleus' at the distances Re ^ = 1. 
In the zero-order approximation we can assume that = and Re = I. This leads us 
to the hydrogen-like quasi-atoms ab/ie and negatively charged ions ab/iee. However, even 
in this approximation we cannot ignore correlations between different particles, and first of 
all electron-muon correlations. In the next order approximation when i?^ > the situation 
becomes more complicated. This requires an additional flexibility of the variational anzatz 
which is used to construct trial wave functions. In all five-body calculations performed 
in this study we have used the variational anzatz which is written in the ten-dimensional 
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gaussoid functions (or five-body gaussoids) of the ten relative coordinates Vij 

N (45) 

^L=0 = (l + 5a6ea6Pa6)(l + P45)E E Ck ■ CXpi-a^f (3) 



k=l fe-)=(12) 



where Ck are the hnear coefficients (or hnear variational parameters), while aij\= aji) are 
the optimized non-linear parameters. It is assumed in Eq.dH]) that i < j. The notation 
Vab means the permutation operator of the two particles a (particle 1) and b (particle 2). 
Additional phase Sab^ab = ^^ab, where 6km is the Kronecker delta, is needed to provide the 
correct permutation symmetry in those systems where a = b. The analogous operator V45 is 
a pair permutation of the two identical electrons (particles 4 and 5) in the five-body abfiee 
ion. The overall electron symmetry of the trial wave function, Eq.([3]), corresponds to the 
singlet electron state (s). 

The results of numerical calculations of the ground state energies of the symmetric five- 
body ions fi~ e~ e~ , d'^ fi~ e~ e~ and t^t^ fi^e^e^ are shown in Table I. In these calcu- 
lations we have used 400 basis functions in Eq.([3]), i.e. = 400. The same Table contains 
the total variational energies of the ground states of the symmetric four-body quasi-atoms: 
p'^ p"^ fi~ e~ , fi~ e~ and t^t^fi^e' (see Section II). As follows from Table I the total 
energies of the five-body ions p'^p~^fi~e~e~,d'^d~^fi~e~e~ and t^t^ ^~e~e~ are lower, i.e. 
they are more negative than the total energies of the corresponding four-body quasi-atoms 
p^p'^fi~e~,d'^d'^fi~e~ and t^t^ji^e~. Note again that in our four-body computations we 
have used the trial wave functions with A^ = 400 and 600 terms and all total energies of 
the five-body ions have been determined with the use of A^ = 400 basis functions, Eq.([3]). 
This means the absolute stability of the five-body ions a+a^yU^e^e^ against dissociation as 
represented by Eq.([2]). By performig analogous calculations for the corresponding five- and 
four-body non-symmetric ions a^b^ fi^e~e~ and quasi-atoms a^b^^~e~ one finds that these 
ions are also stable. 

Let us compare our results obtained for the five-body a^b^ fi^e^e^ ions with the known 
total energies of the hydrogen atom H and hydrogen ion H~ with the infinitely heavy nucleus. 
Traditionally in atomic physics such systems are designated as the °°H atom and °°H~ ion, 



respectively. The best-to-date variational energy obtained for the hydrogen °°H ion {4] is 
-0.527 751 016 544 377 196 589 733 a.u., where 21 - 22 decimal digits are stable. The total 
energy of the non-relativistic °°H atom is -0.5 a.u. (exactly). Therefore, the corresponding 
binding energy of the negatively charged hydrogen ion is ~ -0.027751 a.u. Approximately 



the same binding energy is obtained for each of the five-body ions a"'"a+/i~e~e~ considered 
here. In reahty, one finds the two following differences between the hydrogen °°B.~ ion 
and five-body a^h^ ^~e~e~ ion: (1) the mass of the central quasi-nucleus a^b^fi~ is always 
finite, and (2) there are correlations between electrons and heavy particles from the central 
quasi-nucleus. It is clear that the electron-muon correlations play the leading role here. In 
general, the contribution of such correlations slightly decreases the actual binding energy of 
the five-body a~^b^fi^e^e^ ion. 



IV. HYPERFINE STRUCTURE OF THE GROUND STATES OF THE FOUR- 
BODY QUASI-ATOMS 

Now, consider the hyperfine structure splitting of the ground states in the four-body 
quasi-atoms a~^b^fi^e^ and five-body negatively charged ions a~^b^ fi~e~e~ . Note, that the 
ground state of the five-body a+6+yU~e~e~ ion is the electron singlet state. The total electron 
spin of such a state equals zero exactly. Therefore, the electron spin does not contribute 
to the overall hyperfine structure. Briefly, we can say that the hyperfine structure of the 
five-body ion a^b^ ^~e~e~ exactly coincides with the hyperfine structure of the three-body 
central 'quasi-nucleus' a^b^fi^. The hyperfine structures of the ground S{L = 0)— states in 
six muonic molecular ions a~^b~^fi~ have been analyzed recently in [Sj]. 

In contrast with the five-body (two-electron) ions a'^b'^ fi~e~e~ , the corresponding four- 
body quasi-atoms have very interesting hyperfine structure. Below, in this study we restrict 
ourselves to the consideration of the p^d^ ^^e" and p^t^fi~e~ quasi-atoms. The proton, 
triton, muon and electron have spin equal |, while deuteron's spin equals 1. The dimension 
of the total spin-space for the pdfie quasi-atom is 2x3x2x2 = 24. The representation 
of the rotation group acting in this spin space is represented as the direct sum of the 
corresponding irreducible representations. The ranks of such irreducible representations 
essentially determine the hyperfine structure, while the differences between hyperfine energy 
lev.b ca. be obtained by =olvmg the co.e.ponding eigenvalue problem (all details can be 
found in [71]). The Hamiltonian {AH)hs. which represents the spin-spin interactions takes 
the following form (in atomic units) (see, e.g., js]) 

{AH),... = ^a^^(5(r,.))(s, ■ s,) + ^a^f^(5(r,,))(s, ■ s,) 
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+^a2M£(5(r,^))(s, . s,) + ^a2^(5(r,,))(s, ■ s,) (4) 
3 nip 6 

2 

where a = |^ is the fine structure constant, and nip are the muon and proton masses, 
respectively. The factors g^, gp, gd and ge are the corresponding (7— factors. The expression for 
{AH)h s^ is, in fact, an operator in the total spin space which has the dimension {2sp+l){2sd+ 
l)(2s^ + l)(2se + l) = 24. In our calculations we have used the following numerical values for 
the constants and factors from Eq.(jl]): a = 7.297352586 ■ 10~^, nip = 1836.152701me, = 
206.768262me,^e = -2.0023193043718 and g^ = -2.0023218396. The ^-factors for the 
proton and deuteron are deteremined from the formulas: gp = ^ and gd = where Aip = 
2.792847386 and Aid = 0.857438230 are the magnetic moments (in nuclear magnetons) of 
the proton and deuteron, respectively. The spin of the proton and deuteron is designated in 
Eq.(jlD as Ip= I and Id = 1. 

The energy levels of the hyperfine structure and the corresponding structure splittings 
(differences between energy levels) can be found in Table II. Note that these values are 
usually given in MHz, while the values of {AH)h.s. which follow from Eqs.dl])) are expressed 
in atomic units. To re-calculate them from atomic units to MHz the conversion factor 
6.57968392061 -10^ MHz/a.u. was used. The general hyperine structure of the ground 
state of the pdfie quasi-atom follows from Table II. The total spin function has twenty four 
components. It can be represented as a direct sum of the following irreducible components 
(from the top to the bottom): (1) J = | (four states), (2) J = | (six states), (3) J = | 
(two states), (4) J = | (four states), (5) J = | (two states), (6) J = | (four states) and (7) 
J = 1 (two states). The total final dimension of this direct sum is: 4 + 6 + 2 + 4 + 2 + 4 + 
2 = 24 as expected. The hyperfine structure splittings, i.e. the differences A, are shown in 
Table II. There are three large and three small A. In particular, the values of A12, A34 and 
Ag? are small. Very likely these differences are related with the protium-deuterium spin- 
spin interaction which is very small by its absolute value, since the expectation value of the 
proton-deuteron delta-function is very small On the other hand, the values of A23, A45 
and A56 are relativey large. They are related with the electron-nuclear, electron-muon and 
muon-nuclear spin-spin interactions. A slightly more contrast picture can be obtained for 
the ground states in other muonic quasi-atoms, e.g., dt^e and tt/xe. We have determined 
the hyperfine structures for the ground states in all six muonic quasi-atoms mentioned in 
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Pj. The hyperfine structure of the pt/j^e quasi-atom can be found in Table III. In this case 
gt = where Aip = 2.792847386 and It = |. The hyperfine structure of the pdfie,ptfie 
and other similar quasi-atoms which contain Coulomb three-body quasi-nuclei ab/j, must be 
confirmed in future experiments. 

V. CONCLUSION 

Thus, we have illustrated by numerical means the stability of the bound five-body ions 
abfiee in their ground states. Some of the properties of these five-body ions are also pre- 
dicted. The hyperfine structure of these ions is identical to the hyperfine structure of the 
corresponing muonic molecular ions abfi. We also investigate the hyperfine structure of the 
four-body muonic quasi-atoms abfie which is substantially more complicated. To obtain 
this hyperfine structure splitting we have diagonalized the Hamiltonians of the spin-spin 
interactions which include six different terms and overall dimension up to 36 (for the dd/j^e 
and {ddfi)*e quasi-atoms). The hyperfine structure of the pd/ie and ptfie quasi-atoms is 
investigated in details. 
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TABLE I: The total non-relativistic total energies E of the ground S{L = 0)— states in the five-body 
ions aajiee and four-body quasi-atoms aa/xe (in atomic units). 







ppiiee 


ddfiee 


ttfiee 


E{N 


= 400) 


-102.75052695 


-110.34390610 


-113.49923563 






ppfxe 


ddfj,e 


ttjie 


E{N 

E{N 


= 400) 
= 600) 


-102.72331085 
-102.72332519 


-110.31683802 
-110.31684129 


-113.47266365 
-113.47266986 






pdfxe 


ptjie 


dtfxe 


E{N 
E{N 


= 400) 
= 600) 


-106.51226911 
-106.51233102 


-107.99443693 
-107.99453043 


-111.86419287 
-111.86421475 



TABLE IL The hyperfine structure Ei and hyperfine structure splitting Aj(j_i) of the ground 

S{L = 0)— state of the pd/xe quasi-atom (in MHz). 



Ei MHz MHz 

€. 3 1.2521503100-10'^ 

€t_6 1.2517912194-10^ 3.59090572-103 
ej_i 9.3086776002- 10^ 3.209234594- 10^ 
ej_3 9.3043876675- 10^ 4.28993274-10=^ 
ej_i -2.1222389818-10^ 3.052677749-10^ 
e,_3 -2.3095865825-10^ 1.873476007-10^ 
ej_i -2.3100074249-10^ 4.208423455-10^ 

-'~2 
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TABLE III: The hyperfine structure Ei and hyperfine structure splitting Aj(j_]^) of the ground 
S{L = 0)— state of the four-body pt/ie quasi-atom (in MHz). 



Ei MHz Aj(j_i) MHz 

ej=i 1.5082272852-10^ 

ej=2 1.5078349188-10^ 3.9236642227- 10^ 
ej=o 1.6857006656- 10^ 1.3392648522-10^ 
ej=i 1.6776620651-10^ 8.0386005145-10=^ 
€j=i -3.1838284105-10'^ 3.3515946170-10^ 
ej=o -3.1842399042-10^ 4.1149362917-10^ 
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